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Abstract 

Q We present an analysis of time-delayed feedback control used to stabilize an unstable steady 

state of a neutral delay differential equation. Stability of the controlled system is addressed 
by studying the eigenvalue spectrum of a corresponding characteristic equation with two time 
delays. An analytic expression for the stabilizing control strength is derived in terms of original 
system parameters and the time delay of the control. Theoretical and numerical results show that 
the interplay between the control strength and two time delays provides a number of regions 
I in the parameter space where the time-delayed feedback control can successfully stabilize an 

£>. otherwise unstable steady state. 

C\ 1 Introduction 
in 

Control theory plays an important role in engineering, physics and applied mathematics. It is 
a well known fact in dynamical systems theory that any systems can be controlled by applying 
i—l small perturbations to the available so-called control parameter. Using a control theory approach 

it was shown that it is possible to control unstable or even chaotic oscillations in many systems 
[TJ O [2]. This possibility of control is of a paramount importance in many real-life applications, 
where unwanted oscillations can lead to poor performance and even dangerous operating regimes 
[H [5] . An excellent overview of recent advances in the control of nonlinear dynamical systems can 
be found in Ref . (6] . One of possible approaches to control unstable steady states is an adaptive 
stabilization based on a low-pass filter. Details of such a scheme together with examples, including 
stabilization of a driven pendulum, Lorenz system and an electrochemical oscillator are given in 
Refs. |S]. Another control technique, which is quite often used in applications (e.g. optical 
experiments), and may even be simpler to realize than a low-pass filter, is time-delayed feedback 
control (TDFC), or time-delay autosynchronization first introduced by Pyragas in [9|; a survey of 
the developments of this control method is given in Ref. [TU]. This type of control employs the 
difference between the state of the system at the current time and the state of the system some 
time r ago. In the problem of stabilizing periodic orbits, the time delay r is usually chosen to be 
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the period of the unstable periodic orbit (UPO). This makes the TDFC non-invasive as it vanishes 
as soon as the target state is reached, and it does not require a lot of information about the system 
in general. TDFC has been successfully applied to many physical problems such as, for example, 
controlling unstable motion in electronic chaotic oscillators fT2"[ [TB"] , control of spatio-temporal 
UPOs in coupled reaction-diffusion models [TH [151 US E] • Further modifications of TDFC include 
extended time-delayed feedback control, where multiple time delays are considered in the feedback 
loop [13], unstable feedback controllers |18| . and spatio-temporal filtering [19J. 

Besides control of unstable periodic orbits, TDFC and its multiple time delay extension can be 
successfully used to control unstable steady states [201 ED [221 [23]. In this case, there is an extra 
flexibility in varying the time delay, which in the case of controlling UPOs is required to coincide 
with a period of an UPO to ensure that the control is non-invasive. While TDFC is often unable 
to stabilize a state having a single real positive eigenvalue [21] (albeit this is not generally true 
for UPOs [25]), it is quite successful in stabilizing steady states with a pair of unstable complex 
eigenvalues |22l 123] . Asymptotic analysis of the appropriate characteristic equation shows that, 
provided the steady state is close enough to the instability threshold, the TDFC provides efficient 
stabilization even for relatively large time delays |26j. 

In this paper we use time-delayed feedback control to stabilize an unstable steady state of a 
neutral delay differential equation (NDDE). A delay differential equation is called neutral if it 
contains a time delay in the highest derivative involved. This type of equations arises in numerous 
physical and engineering application, for example, hybrid testing |271 128] . chaotic oscillations in 
transmission lines |29j and torsional waves of a drill string [30J. Due to the original time delay 
present in the system, its steady states may become unstable through a Hopf bifurcation, and 
adding a time-delayed feedback control will stabilize the system again. The interplay between the 
control strength, the internal and the external time delays, leads to the existence of various regimes 
in the parameter space where the steady state of the controlled system is stable. At the same time, 
due to essential instability in the spectrum, it is impossible to achieve controllability when the value 
of the coefficient multiplying the highest order time-delayed derivative in the equation exceeds one. 

The paper is organized as follows. We first present a stability analysis of an NDDE with a 
single time delay. Then we introduce time-delayed feedback control to stabilize an unstable steady 
state of the original system. We show that by varying the control strength and the time delay 
it is possible to expand the stability region for a range of time-delays. To better understand the 
influence of TDFC on the stability of the system, we numerically compute the eigenvalue spectrum 
of the fixed point. The paper concludes with a summary. 

2 Stability analysis 

Consider the following NDDE derived in Kyrychko et al. [27], which describes a pendulum-mass- 
spring-damper system: 

z(t) + 2(z(t) + z{t) + pz(t - r) = 0, (1) 

where dot means differentiation with respect to time t, r is the time delay, £ and p are model 
parameters corresponding to the damping constant and the ratio of the mass of the pendulum to 
the mass of the mass-spring-damper, respectively. We will use this equation as a prototype model 
for the linearization of other neutral DDEs (see [29] and [30J for examples) near their steady states. 
This equation has a single steady state z* = 0, whose stability is determined by the roots of the 
corresponding characteristic equation 

A 2 + 2CA + 1 + p\ 2 e~ Xr = 0. (2) 
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Figure 1: Stability boundary of the system ([!]) for £ = 0.1. The steady state is stable in the shaded 



area. 



In the case when the parameter p in Eq. ([!]) exceeds unity, the steady state is unstable for any 
positive time delay r. As it was shown in Kyrychko et al. [27], for \p\ < 1 the steady state 
switches stability as the time delay is varied. The stability boundary can be parametrized by the 
Hopf frequency ui. In this case, the critical time delay and the critical value of p at the stability 
boundary are given by 



UJ 
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1) 2 + 4C 2 ^ 2 , 



where k = 0, 1, 2, ... and Arctan denotes the principal value of arctan. Figure [T] illustrates the 
dependence of the critical value of p on the time delay r which ensures the stability of the steady 
state. It is noteworthy that if \p\ > 1, the steady state is unstable for any positive time delay r 
due to the so-called essential instability [3T] . On the other hand, if \p\ < 1 and £ > l/\/2, then the 
steady state is asymptotically stable for any positive time delay r. For ( < l/v2, there is a lower 
bound on the value of p m i n = 2£y / l — £ 2 , so that asymptotic stability is guaranteed for all r > 
provided p < p m ; n [27] . 

In the area above the boundary curve in Fig. [TJ the steady state of Eq. is unstable, with 
a pair of complex conjugate eigenvalues in the right complex half-plane. The aim of this study is 
to stabilize this unstable steady state by means of time-delayed feedback control. To this end, we 
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Figure 2: (Color online) (a) Stability boundary of the controlled system for t\ = r 2 = r, £ = 0.1 
and different control strengths: K = (solid), K = 0.1 (dashed) and K = 0.2 (dotted). The steady 
state is stable below the curves, (b) Minimum stabilizing control strength K as a function of p and 
ti = r 2 = r. 



replace the original system ([I]) by its modified version 

z(t) + 2(i(t) + z(t) + - n) = #[z(i) - z(i - r 2 )], 



(3) 



where if > is the strength of the control force, and r 2 > is the time delay of the control. This 
modified system has the same steady state as the original system for any time delay r 2 , but now 
the stability of this equilibrium is determined by the roots of the modified characteristic equation: 



A z + 2CA + 1 + pX z e- XTl =K(1 



-AT2 



(4) 



If K = or r 2 = 0, we recover the original uncontrolled model with an unstable steady state. The 
goal is now to find a relation between K and t 2 in terms of the original system parameters p, 
and ti, so that in the absence of control there is a pair of unstable complex conjugate eigenvalues, 
while in the presence of control all the eigenvalues are stable. 

First we consider the case when the time delay of the control coincides with the time delay of 
the original system r 2 = t\ = r. Let us determine the stability boundary Re(A) = 0. Looking for 
solutions of Eq. Q in the form A = iuj and separating real and imaginary parts, we find 



-oj 2 + 1 
2Qu = - 



- K = (pu 2 — cos cjt, 
{poj 2 — if) sin ujt. 



(5) 



Squaring and adding the two equations in the last system gives 
The last equation can be solved for p as follows: 



P 
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(-6J 2 + 1 - K) 2 + 4( 2 u 2 + K 
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Also, by dividing the second equation of the system ([5j by the first, the critical time delay can be 
expressed as 

1 



to 



2(lu 

Arctan — = ± nk 

uj z — 1 + K 



k = 0,1,2, 



Now, for each fixed value of the control strength K and the damping £, we have a stability boundary 
in the (r,p) plane parametrized by a Hopf frequency uj. Figure [2]^a) shows that by increasing the 
control strength K, it is possible to raise the stability boundary, and thus improve overall stability 
of the system. This also means that those points which were unstable in the uncontrolled system, 
are now stable as they are moved by the control action to the area below the new stability boundary. 
It is worth mentioning that if K < 0, this actually lowers the stability boundary, and hence reduces 
the stability. 

So far, we fixed the control strength K and studied the stability boundary as a curve in the (r,p) 
plane. At the same time, one can consider p and r fixed (this corresponds to a single unstable fixed 
point of the uncontrolled system), and then solve the system ^ for a tuple (K,oS). This would 
give a full picture of controllability of the system by providing the value of the minimum control 
strength K required to stabilize the steady state for the given values of p and r. The results of this 
computation are shown in Fig. [2^b), where the steady state is stable above the surface. This plot 
suggests that the more unstable is a steady state (i.e. the higher it is above its stability boundary 
in (r,p) plane), the higher should be the control strength K required to stabilize this steady state. 

When T2 7^ ti, the characteristic equation Q has two time delays simultaneously present, which 
significantly complicates a stability analysis. Several approaches have been recently put forward 
to study stability of equations with multiple time delays (20J EE]- Beretta and Kuang [3l] have 
developed geometric stability switch criteria for equations with delay-dependent parameters, and 
this method can also be used to analyze systems with two time delays. Sipahi and Olgac |32j 
have suggested a systematic procedure of finding eigenvalues of characteristic equations for systems 
with multiple time delays by means of a substitution, which replaces an original transcendental 
characteristic equation with a polynomial. For the purposes of this paper, it is convenient to use 
a parametrization of the stability boundary in the parameter plane of two time delays, as outlined 
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Figure 4: (Color online) Real part of the leading eigenvalue of the characteristic equation Q versus 
delay times t\ and T2 for p = 0.4 and £ = 0.1. (a) K = 0.1. (b) K = 0.3. The color code denotes 
the value of Re(A); only those parts are shown where it is negative. The colored areas in both plots 
correspond to the regions where the steady state is stable. 



in Gu et al. [33] . The idea is to rewrite the characteristic equation Q in an equivalent form 

1 + ax (A)e~ Ari + a 2 e- XT2 = 0, (6) 

where 



oi(A) 



1 - K + 2CA + A 2 ' 



02(A) 



K 



I- K + 2CX + X 2 ' 



The stability border can now be parameterized by a Hopf frequency uj £ f2, where Q = IJi=i 
consists of a finite number of intervals of finite length. The critical time delays at the stability 
boundary in the (ti,7^) plane are then given by 



n 



Arg[ai(iw)] + (2u - 1)tt ± 61 



u = uf ,Uq + 1,Uq + 2, ... 



T2 



Arg[a 2 (ia;)] + (2u - 1)tt ± 0i 



(7) 



where 61,62 £ [0, 7r] are calculated as 



6\ = arccos 



12 = arccos 



1 + |ai(iw)| 2 - \a 2 {iu)\ 2 
2|ai(ia;)| 

1 + |a 2 (iw)| 2 - \ai(iui)\ 2 
2|a2(«w)| 



and Uq and Vq are determined as the smallest positive integers such that the corresponding 



r"° , Tg , T2 are all non-negative. 
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Figure 5: (Color online) (a) Real part of the leading eigenvalue of the characteristic equation Q 
for p = 0.4, C = 0.1, and t\ = 5. b) Dependence of the real part of the leading eigenvalue on the 
control time delay T2 for K = 0.15 (solid), K = 0.28 (dashed), and K = 0.36 (dotted). 



Figure [3] shows the stability boundaries of the controlled system in the parameter space of the 
time delays T\ and T2, as given by the formulas ([7]). We note that when K = 0.1, the stability 
boundaries are represented by zig-zag curves with vertical axes. For larger values of the control 
strength K, the stability boundaries become closed curves. The steady state is stable inside the 
shaded areas. 

The control strength K appears as a parameter in a\ and a>2, and thus determines the values of 
T\ and T2 at the stability threshold. At the same time, in order to determine the controllability of 
the system it is more convenient to fix particular values of the parameters t\ , p, and £ of the original 
system, and then study the stability border in the (t2,K) space. Figure [5^a) shows how the real 
part of the leading characteristic eigenvalue depends on K and T2, and in Fig. [5^b) for illustrative 
purposes we present three possible regimes for different values of the control strength K. One can 
observe that while for small K the steady state is unstable, it becomes stable starting with some 
K, and the larger K is, the more negative the real part of the leading eigenvalue becomes. It is 
clear from Fig. [5^a) that the stability islands in (t2, K) parameter plane are finite, and hence if the 
control strength is very large, the steady state will be unstable again. 

If the parameters are taken closer to the stability boundary of the uncontrolled system, a much 
smaller control strength K is required to stabilize an unstable equilibrium. 

To deepen the understanding of the stability changes depending on the relation between the 
two time delays, we have computed the maximum of the real part of the characteristic eigenvalues, 
which is shown in Fig. [4] wherever it is negative. To compute the real part of the leading eigenvalue 
of the characteristic equation Q we have used pseudospectral differentiation methods [351 GS EZ] • 
Note some of the very small stable domains (with Re(A) very close to zero) shown in Fig. [3] are not 
fully resolved numerically in Fig. H| In this case all stability boundaries are closed curves shown in 
Fig.[6]for fixed K = 0.1 and K = 0.3, and the steady state is stable in the shaded area outside those 
curves. The corresponding plot of the real part of the leading eigenvalue, shown in Fig. [7j reveals 
the existence of several regions of values of T2, for which the steady state can be stabilized for the 
same value of the control strength K and the original time delay t%. If we fix n and consider the 
plane of control parameters K and T2 , the dependence of the leading eigenvalue of the characteristic 
equation on these parameters is qualitatively the same as the one for p = 0.4, whereas the stability 
regions are more pronounced, see Fig. [8j In comparison to the case of p = 0.4, there is more than 



7 





Figure 6: Stability boundaries of the characteristic equation Q for p = 0.25 and £ = 0.1. a) 
K = 0.1, b) K = 0.3. Shaded area corresponds to stable steady state. 

one interval of time delay T2 , for which stabilization of the steady state is achieved for a fixed value 
of the control strength K. Note that similar behavior as in Figs. [6] and [7] was found for an electronic 
oscillator where two time delays were added in the form of two Pyragas-type time-delayed feedback 
control terms |21j . 

3 Conclusions 



For many physical, biological and engineering problems an important practical challenge is the 
control of stability of the steady states. While several techniques are readily available, not all of 
them are equally robust and can be implemented experimentally. In this respect, time-delayed 
feedback control has been demonstrated to be an efficient tool for stabilization of unstable periodic 
orbits and unstable steady states. 

The instability of the steady states occurs quite often in time-delayed systems, where the time 
delay can destabilize the steady state through a Hopf bifurcation. In this paper, we have shown how 
the time-delay feedback approach can be used to stabilize unstable steady states in a neutral delay 
differential equation. For a given value of the control strength, it is possible to find an analytical 
parametrization of (in)stability boundaries in terms of two time delays: t\ (time delay of the original 
system) and ti (time delay of the control force). Quite naturally, the more unstable is a steady 
state, as determined by the real part of the leading eigenvalue of the corresponding characteristic 
equation, the higher should be the control strength required to stabilize it for a fixed value of the 
control time delay T2- At the same time, if the steady state is close to a Hopf bifurcation, it is 
easier to control it. In the latter case there is a larger number of islands in a parameter space where 
stabilization is successfully achieved. 

It is noteworthy that the time-delayed feedback control scheme fails to achieve controllability 
of the steady state if \p\ > 1 independently on the intrinsic time delay and control parameters. 
This further highlights the difference between stabilizing unstable steady states in neutral and 
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Figure 7: (Color online) Real part of the leading eigenvalue of the characteristic equation Q versus 
delay times t\ and T2 for p = 0.25 and C, = 0.1. (a) X = 0.1. (b) K = 0.3. The color code denotes 
the value of Re(A); only those parts are shown where it is negative. The colored areas in both plots 
correspond to the regions where the steady state is stable. 

non- neutral time-delayed systems. 
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